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Abstract 

This note answers one question in [3], concerning the connected allo- 
cation for the Poisson process in R 2 . The proposed solution makes use 
of the Riemann map from the plane minus the minimal spanning forest 
of the Poisson point process to the halfplane. A picture of numerically 
simulated example is included. 

1 The problem 

Let X C M d be a discrete set. We call the elements of X centers, the elements 
of M. d - sites, and we write C for the Lebesgue measure in K 2 . An allocation of 
M. d to X with appetite a £ [0, oo] is a measurable function -0 : M. — > X U {oo, A} 
such that A^HO < a for all £ E X and £V -1 (A) = 0. We call V _1 (0 
the territory of the center £. A center £ is sated if £i/> _1 (£) = a and unsated 
otherwise. A site x is claimed if VK^) £ an d unclaimed if VK 2 -) = 00 • The 
allocation is undefined at a; if ^p{x) = A. 

One particular question we're interested in is the following [3]: Is there a 
translation- equivariant allocation for the Poisson process of unit intensity in 
the critical two-dimensional case (d = 2 and a = 1), in which every territory 
is connected? The goal of the present note is to describe one such allocation. 
Note: for a far more general approach to this question in d > 3 see [5] . 

2 Construction of a connected allocation 

Let X be a realization of Poisson process (of intensity 1) in R 2 . 

Let the minimal spanning forest T be the union of edges e = (x, y), x, y € X, 
such that there is no path from x to y in the complete graph on X with all edges 
strictly shorter than e (see [J] for other definitions and references on the subject). 

For a 2-dimensional Poisson process Alexander [T] proved that T is a.s. a tree 
with one topological end (a topological end in a graph is a class of equivalence 
of semi-infinite paths modulo finite symmetric difference). Thus the domain 
D = M 2 \T is connected, and by the Riemann mapping theorem can be mapped 
conforamlly to the upper half-plane EI = {z 6 C : 3(z) > 0}. Moreover, we can 
choose such a mapping / which also sends infinity to infinity. (Note that any 
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two such functions differ by a conformal automorphism of H which invariates 
infinity, i.e. by an affinc mapping z — > az + b, a, b <G M, a^O). 

The mapping / cannot be extended unambiguously to include T. Instead, 
consider the inverse mapping g — f" 1 from the open halfplane to D, then extend 
g to the real axis by continuity, obtaining a surjection ~g : H — > M 2 . Some points 
of X may have more than one preimage under g, more exactly, a point has as 
many preimages as it's degree in T. For those points we'll have to split their 
appetite between these preimages in some way, e.g. proportionally to the angle 
of the corresponding sector. Let (x' n , a' n ) be the resulting points and appetites. 
In the following we call the countable set {(x' n ,a' n ),n G Z} the image of X 
under f. 
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Figure 1: The image of X under /; the unit appetite of the center A is split 
into three unequal parts. 

Now consider the following stable allocation procedure (see Appendix II for a 
formal definition): 

• Each center starts growing a ball centered in it. All the balls grow simul- 
taneously, at the same linear speed. 

• Each center claims the sites captured by it's ball, unless the site was 
claimed earlier by some other center. 

• Once the center becomes sated (the measure of it's territory reaches it's 
appetite), the ball stops growing. 

Note that from our choice of function / it follows that the set {x' n ,n £ Z} is 
locally finite (i.e. has no accumulation point other than infinity), so the stable 
allocation procedure is well defined. 

Applying this procedure to the image of X under / in the half-plane H, with 
respect to the Euclcdian metric in H and the image A of the Lebesgue measure 
C under / yields certain allocation tpn :B^IU {oo, A}. 

Lemma 2.1 The allocation ipu thus constructed has connected territories and 
is invariant under affine transformations o/H. 

Proof. The affine invariance follows immediately by construction. 

To show that each center has connected territory, consider the following. Let 
A be a center, necessarily located at the boundary of H, and let Da = ^^(A) 
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be it's territory. Consider the ray AA' , perpendicular to R. If some point x on 
A A' doesn't belong to Da, that's because the center got sated before reaching 
x, so no further point on AA' belongs to Da- Thus the intersection of Da with 
A is a segment AK. 
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Figure 2: Every point of Da can be reached from A 

Now consider an arc of a circle centered in A and intersecting AK in some 
point Y, and follow this arc to the right (or left) from Y . Once we meet a point 
Z Da, there must be a center B such that BZ < AZ, so no further points on 
the arc belong to Da- Thus every point of Da can be reached from A, and Da 
is connected. 

Now apply the inverse map / _1 to ipH- i-e- let ip = f~ 1 ipnf- Clearly, ip 
is an allocation of R 2 to X with appetite a, and from the previous lemma it 
doesn't depend on the choice of /. 

Lemma 2.2 The allocation ip is translation-invariant, i.e. if t : R 2 — > R 2 is 
a translation, and ip' it the allocation or R 2 to tX constructed using the above 
procedure, then ip' — T~ x ipr . 

Proof. Let X' = rX. Clearly, the minimal spanning forest is translation- 
invariant, i.e. V := MSF(X') = t ■ MSF(X). Now /' = / r^ 1 is a conformal 
function that maps R 2 \T' to H, and the image of X' under /' coincides with 
the image of X under /, thus ip' = f'~ 1 ipHf' = T~ l f^ipiifT — t _1 ^t. 

Lemma 2.3 Under the allocation ip every center is sated a.s. 

Proof. First, it follows from the Galey-Shapley algorithm (see Appendix II 
below), that no stable allocation may have both unclaimed sites and unsated 
centers. By ergodicity, the existence of unclaimed sites is a 0/1 event; thus we 
may assume that ipu, and therefore ip, have no unclaimed sites. 

Now Lemma 16 in [5] states that for any translation-invariant allocation ip 
and any r > 

P{\iP(0)\ <r} = E* C(ip- x (0) n B(0, r)) , 

where E* denotes the expectation with respect to the Poisson process condi- 
tioned to have a center in 0. Taking r — > oo yields 

P{0 is claimed} = E* Cip- l {Q). 
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But wc assumed the probability on the left to be equal to 1, and £^i> -1 (0) < 1 by 
construction (a center cannot allocate more than it's appetite), thus the lemma 
follows. 

We summarize the above lemmas in the following 

Theorem 1 The allocation ip thus constructed has connected territories, is 
translation-invariant, and every center is sated a.s. 

3 Further questions 

Geometry of territories. Is it true that the every territory V^ 1 ^) is bounded 
a.s.? The territories of tp are connected, but, a priori, may be not simply 
connected (more exactly, their closures may not be simply connected). Can the 
construction of tp be modified to assure this? 

Point process on R. What can be told about the image of X under the 
Ricmann map, as a point process on 7?? Simulations suggest that this process 
should have extremely high variation of point density; even for small polygons 
of 30-40 points, the distances between the images of consecutive points may 
vary from rj 1 to ~ 10~ 16 . 

Fixing the Riemann map. The Riemann map / is unique up to a conformal 
automorphism of H , z — > az + b. Is it possible to pick (deterministically) one 
particular mapping from this class? If omitting the translation-invariance, one 
obvious way to do this is to pick a mapping / that sends to i. If doing this in 
translation-invariant manner, there is obviously no way to fix the shift b (since 
this would imply the deterministic translation-invariant choice of single point 
from the Poisson process). But does there exist a way to pick the scale, i.e. fix 
the parameter a? 
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Appendix I: An obligatory picture 




Figure 3: An approximation of allocation for 164 points in the unit circle 

On the above picture we took 164 uniformly distributed points inside the 
unit circle, and mapped each component Kj of their convex hull minus their 
MST to a region of H, so that all the internal edges are mapped to M. Then for 
each component we allocated it's whole area to the bounding corners, choosing 
the appetites proportional to the angular measure. 

The mapping was approximated with a version of the geodesic algorithm, 
adapted for domains with non-Jordan boundary (or more exactly for domains 
with boundary consisting of a Jordan curve plus a number of internal trees); for 
discussion of numeric algorithms for quasi-conformal mappings see e.g. [5] and 
references therein. 

Since the ratio area/{sum of appetites} varies in different components Kj, 
not all domains have exactly the same area; there are also some artifacts related 
to the limited numeric accuracy. Nevertheless, we expect the picture close to 
the center to be an accurate approximation of the construction described in this 
note. 
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Appendix II: site-optimal Gale-Shapley algorithm 



(The description of the G-S algorithm below is taken from J3jj except for the last 
step, when we have to be more careful with the definition of if) onW in order to 
assure the connectedness of the territories.) 

Let W be the set of all sites, equidistant from one or more centers. Since 
the set of centers is countable, W has Lebesgue measure null. 

We construct tp on R 2 \M^ by means of a sequence of stages, where stage n, 
n = 1, 2, 3, . . . consists of two steps: 

a) Each site x ^ W applies to the closest center, which has not rejected x at 
any earlier stage. 

b) For each center £, let A n {£) be the set of sites which applied to £ on step 
a of stage n, and define the rejection radius as 

r n (0 = inf{r : C(A n ^) n r)) > a}, 

where the infimum over the empty set is taken to be oo. Then £ shortlists 
all sites in A n (£) fl B(£,r„(£)), and rejects all sites in A n (£)\B(£,r n (£)). 

Now either x is rejected by every center (in order of increasing distance 
from x), or i is shortlisted by £ for some center £ at some stage n. In the 
former case we put ip(x) = oo (so x is unclaimed), in the latter case we put 

$(x) = e 

Finally, for x £ W put ip(x) = £, if £ is the only center such that x £ dip 1 (£), 
and put ip{x) — A (undefined) otherwise. 
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